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I .  INTRODUCTION 


The  mechanism  of  solid  propellant  combustion  is  commonly  idealized  as  a 
single  exothermic  pyrolysis  reaction  occurring  at  the  surface  followed  by  a 
single  exothermic  chemical  reaction  in  the  gas  phase.  Additional  assumptions, 
such  as  constant  thermal  conductivity,  and  unit  Lewis  number,  are  also  usual. 

The  appropriate  equations  and  boundary  conditions  are  derived  elsewhere.* 

A  comparison  of  the  solution  to  these  equations  with  experimentally  measured 
burning  rates  should  provide  information  as  to  the  correctness  of  the  assumed 
mechanism.  In  practice  one  of  the  difficulties  with  this  line  of  reasoning 
(others  are  iiscussed  in  Ref.  2)  is  that  these  equations  do  not  have  analytic 
solutions  and  indeed  can  sometimes  be  difficult  to  solve  numerically.  Generally 
this  problem  has  been  dealt  with  by  making  further  simplifications  in  the 
equations  to  obtain  a  convenient  computational  algorithm.  The  object  of  this 
paper  is  to  test  the  accuracy  of  these  algorithms  by  direct  comparison  with 
numerical  solutions. 

It  has  previously  been  shown1  that  most  models  addressing  the  stated 
mechanism  can  be  catalogued  by  the  assumed  manner  of  heat  release  in  the  gas- 
phase.  These  methods  of  specifying  the  heat  release  are  generally  appropriate 
to  complementary  ranges  in  the  kinetics  parameters.  Since  no  general  agreement 
exists  on  the  values  of  these  parameters  for  any  propellant,  we  have  developed 
several  data  sets  which  cover  a  sufficient  range  of  variation  to  probe  the 
strengths  and  weaknesses  of  all  the  model  algorithms.  The  reasoning  behind 
these  data  sets  is  given  in  a  companion  paper2  but  the  values  are  summarized 
here  in  Table  1  for  convenience. 

II.  DESCRIPTION  OF  PROBLEM  AND  NUMERICAL  SOLUTION 
A.  Basic  Idealization. 

The  mechanism  used  in  this  study  is  common  to  more  models  than  probably 
any  other.  Along  with  the  attendant  conservation  equations,  it  is  thoroughly 
discussed  in  Reference  1.  In  brief,  it  is  assumed  that  propellant  molecules 
"A"  pyrolyze  at  the  solid  surface  into  gaseous  products  "B"  and  "C"  at  the 
mass  burning  rate  M  (gm/cm2s)  which  depends  on  the  surface  temperature  Tg 
according  to 

M  =  M  exp(-E  /RT  )  .  (1) 

o  s  s 

(Symbols  are  identified  in  the  List  of  Symbols  and  in  Reference  1).  This 
surface  reaction  is  exothermic  by  an  amount  Qs(cal/gm).  The  pyrolysis  product 
B  is  the  reactant  for  the  gas-phase  reaction  of  order  v.  B  is  converted  by 
this  reaction  into  non-reactive  products  C  with  the  release  of  calories 
per  gram  of  B. 

1 Miller ,  M.S. ,  "In  Search  of  an  Idealized  Model  of  Homogeneous  Solid  Propellant 
Combustion ,"  Combustion  and  Flame,  Vol.  46,  pp.  51-73,  1982. 

2 Miller ,  M.S.  and  Coffee,  T.P. ,  "A  Fresh  Look  at  the  Classical  Approach  to 
Solid  Propellant  Combustion  Modeling ,"  to  be  published  in  Combustion  and 
Flame. 
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In  the  pyrolysis  reaction  for  every  unit  mass  of  A  the  fraction  mB  is 
converted  to  B  and  (l-mB-°)  is  converted  directly  to  C.  At  some  temperature 
T,in  the  gas, the  rate  at  which  B  is  converted  to  C  is  given  by 

B  =  (n>BP)V  Ag  expC-Eg/RT)  (2) 

where  mB  is  the  local  mass  fraction  of  B  and  p  is  the  local  mass  density  of 
the  gas:  For  every  mole  of  B  which  reacts  ^  moles  of  C  are  produced. 

B.  Numerical  Solution. 

The  linear  burning  rate  r  is  to  be  calculated  for  a  semi-infinite  block 
of  propellant  whose  surface  is  situated  normal  to  the  x  axis.  Because  the 
propellant  is  homogeneous,  the  burning  is  one-dimensional.  Two  levels  of 
approximation  to  the  conservation  equations  are  discussed  in  Reference  1. 

The  more  general  of  the  two  is  discussed  here  (termed  Level  I  in  Ref.  1). 

The  method  used  to  obtain  a  burning  rate  in  this  paper  has  been  discussed 
schematically  in  Reference  1.  First  the  conductive  heat  feedback,  <j>G,  from 
the  gas  phase  to  the  surface  of  the  solid  is  computed  numerically  for  many 
values  of  Ts  at  a  given  pressure.  A  value  for  Ts  is  then  sought  for  which 
(j)G  is  equal  to  the  heat  flux,  $s,  required  by  the  solid  tOjmaintain  a  mass 
regression  rate  M  when  the  inital  temperature  is  T  ,  i.e., 

♦,  ■  •vvwv-  (3> 

By  using  this  method  the  numerical  functions  4>r(T  ,P)  do  not  depend  on  either 

A  m  vJ  3 

Qs  V 

The  gas-phase  equations  are  solved  by  a  computer  code  designed  to  handle 
premixed,  laminar,  one-dimensional  flames  involving  elementary  chemistry.  » 
This  code  is  based  on  the  package  PDECOL  developed  by  Madsen  and  Sincovec4 5  to 
solve  sets  of  second  order  partial  differential  equations  using  finite  element 
techniques. 

Although  the  equation  governing  mg  can  be  eliminated  using  the  unit  Lewis 
number  assumption,  both  the  mg  and  the  T  equations  are  integrated  by  the  code. 
The  self-consistency  of  the  two  solutions  can  then  be  checked,  as  will  be 
explained  (Eq.  4). 

3 Coffee,  T.P.  and  Heimerl,  J.M. ,  "A  Method  for  Computing  the  Flame  Speed  for 
a  Laminar,  Premixed,  One  Dimensional  Flame,"  Ballistic  Research  Laboratory 
Technical  Report  ARBRL-TR-02212 ,  January  1980.  (AD  A082803) 

4 Coffee,  T.P. ,  "A  Computer  Code  for  the  Solution  of  the  Equations  Governing  a 
Laminar,  Premixed,  One  Dimensional  Flame,"  Ballistic  Research  Laboratory 
Memorandum  Report  No.  ARBRL-MR-03165,  April  1982.  (AD  A114041) 

5 Madsen,  N.K.  and  Sincovec,  R.F.,  "PDECOL,  General  Collocation  Software  for 
Partial  Differential  Equations,  1  ACM  TOMS,  Vol.  S,  pp.  326-351 ,  1979. 
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The  basic  numerical  procedure  is  to  work  with  the  time  dependent  equations 
using  arbitrary  startirg  profiles  for  mg  and  T.  The  cold  boundary  condition  is 
just  T*T  .  For  the  mass  fraction  equation  it  is1 


+o 


m. 


m 


B 


which  expresses  the  required  continuity  of  B  mass  flux  across  the  boundary. 

The  first  derivatives  of  mg  and  T  are  taken  to  be  zero  at  the  hot  boundary. 
Numerically  the  hot  boundary  was  chosen  sufficiently  distant  from  the  surface 
that  the  mg  and  T  profiles  were  essentially  unaffected  by  its  location.  These 
boundary  conditions  were  previously  used  to  model  burner  stabilized  flames. 

The  equations  are  integrated  in  time  until  the  steady  state  solution  is 
reached.  Since  a  finite  element  code  is  used,  the  solution  profiles  are 
represented  by  piecewise  polynomials  and  the  temperature  gradient  at  the  sur¬ 
face  is  available.  As  a  check  the  heat  feedback  is  also  calculated  by  the 
formula 

+G  =  V,(mB°S0)  t4) 


which  follows  from  the  unit  Lewis  number  assumption. 

The  steady  state  solution  is  usually  quickly  reached,  taking  on  the  order 
of  a  second  on  a  CDC  7600.  The  exception  is  when  M  becomes  so  large  (due  to 
high  Ts)  that  the  flame  is  blown  far  from  the  surface.  This  strong  convective 
flow  allows  little  heat  to  be  conducted  back  to  the  solid.  As  a  result  the 
flame  is  barely  stabilized  and  tends  to  oscillate  before  finally  approaching 
the  steady  state. 

As  a  partial  check  on  the  code,  it  was  run  on  a  test  data  set  (Tl)  for 
which  an  exact  analytic  solution  exists  (Eq.  31  of  Ref.  1).  The  numerical 
solution  proved  accurate  to  within  a  small  fraction  of  a  percent  except  when  Tg 
became  large  enough  so  that  the  heat  feedback  was  negligibly  small.  Even  5 
there  the  numerical  solution  was  good  to  a  few  percent  (and  could  be  improved 
by  using  a  finer  grid) .  In  the  subsequent  figures  the  numerical  computations 
will  be  denoted  by  N. 

C.  Analytic  Models  Used  for  Comparison  . 

Most  of  the  models  which  deal  with  our  assumed  mechanism  have  been  recently 
reviewed  and  classified  into  one  of  three  idealizations  of  the  gas-phase  heat 
release.1  The  constant  temperature  reaction  rate  (CTRR)  approach  should  apply 
in  the  limit  of  very  low  Eq.  The  qua si -const ant  heat  release  (QCHR)  model1 
might  pertain  to  moderate  values  of  Eg  whereas  the  delta  function  heat  release 
(DFHR)  models  would  be  appropriate  to  very  high  Eg-  We  shall  in  this  paper 
compare  our  numerical  calculations  to  the  following  specific  models  in  each 
of  these  classes. 

A  solution  to  $g  under  the  CTRR  assumption  is  not  available  for  v=2  in  a 
Level  I  (diffusion  not  neglected)  treatment1;  therefore.  Equation  (33)  of 
Reference  1  is  used  to  compute  $g  for  the  NCI  data  set.  (Eq.  31  of  Ref.  1 
is  used  for  the  v=l  cases).  A  comparis  of  the  '  lutions  with  and  without 
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diffusion  for  the  v=l  cases  shows  that  the  burning  rate  is  not  greatly  affect¬ 
ed  by  this  assumption  so  long  as  the  solution  value  of  Ts  lies  to  he  right 
of  the  peak  in  the  ♦q(T-)  curve.  For  example.  Fig.  5  suggests  that  the  use 
of  the  no-diffusion  4»G  is  probably  justified  for  Tq  greater  than  about  3C0°K. 

It  became  apparent  during  the  course  of  this  work  that,  while  the  QCHR 
model  $g  gives  good  results  for  $g  at  high  Ts  for  some  cases,  at  sufficiently 
low  Ts  the  model  aG  exceeds  the  maxircum  possible  heat  feedback  (Qq^b-0)  . 

This  occurs  because  at  low  Ts  the  characteristic  transport  distance  (X/MC  ) 
becomes  so  large  (as  M  diminishes)  that  the  integrated  heat  release  " 

exceeds  QGmg-o  .  This  problem  is  easily  corrected  by  assuming  a  constant  heat 
release,  q  ,  only  for  x<x^  where  Xj  is  chosen  by  the  energy  conservation 
requirement 


-o  p 

aB  Vo 


qo  dx 


Thus  we  have  the  "constant  heat  release"  (CHR)  model  given  by 
Xq 

♦g  "  m <T  [i-eaqpC-MC  Xj/A)]  (5) 

P 

which  will  be  used  in  the  present  paper  (cf.  Eq.  (46)  of  Reference  1).  At 
high  T  the  QCHR  and  CHR  results  for  merge.  The  unknown  value  of  mB+° 
in  qQ  is  found  by  requiring  self-consistency  between  Equations  (4)  and  (5). 

In  Reference  1  it  was  argued  that  the  BDP  monopropellant  model^,  although 
a  delta  function  model,  is  likely  to  overestimate  <f>G  in  the  limit  of  high  E^ 
because  of  an  improper  choice  for  the  flame  standoff  distance.  This 
assessment  is  borne  out  in  the  calculations  to  follow.  Fortunately  there  does 
exist  an  accurate  high  EG  theory  which  was  advanced  by  Williams8  and  elaborated 
by  Buckmaster,  et.al.  The  theory  is  based  on  the  asymptotic  analysis  of  9 
Bush  and  Fendell9 10  for  laminar  gas  flames.  In  our  notation  the  burning  rate 
for  general  reaction  orders  is  given  by 

a  Dv+lr  v_  2  (v+1)  _  „ 


2v!XPV  A  RV+1C  VTf2 

(Wb°)V+1 


LRTf' 


exp(-EG/RTf) . 


In  this  paper  we  shall  designate  the  theory  by  A. 

6Beok8tead ,  M.W.,  Dew,  R.L.,  and  Price.  C.E.,  "The  Combustion  of  Solid  Mono-  . 
propellants  and  Composite  Propellants^  Thirteenth  Symposium  ( International ) 
on  Combustion,  The  Combustion  Institute,  pp.  1047-1056,  1971. 

?Hermance ,  C.E.,  "A  Model  of  Composite  Propellant  Combustion  Including  Surface 
Heterogeneity  and  Heat  Generation,”  AIAA  Journal,  Vol.  4,  pp.  1629-1657 ,  1966. 

8 Williams,  F.A.,  "Quasi-Steady  Gas-Flame  Theory  in  Unsteady  Burning  of  a  Homo¬ 
geneous  Solid  Propellant.”  AIAA  Journal,  Vol.  11,  pp.  1328-1330,  1973. 

9 Buckmaster,  J.D.,  Kapila,  A.K. ,  and  Ludford,  G.S.S.,  "Linear  Condensate 
Deflagration  for  Large  Activation  Energy  ”  Acta  Astronautica ,  Vol.  3,  pp.  593- 
614,  1976. 

10Bueh ,  W.B.  and  Fendell,  F.E.,  "Asymptotic  Analysis  of  Laminar  Flame  Propagation 
for  General  Lewis  Humbers,"  Comb.  Sci.  Tech.,  Vol.  1,  pp.  421-428,  1970. 


III.  RESULTS 


I 


i 


i 


i 


i 


i 


In  this  section  we  shall  first  show  how  the  model  heat  release  and 
temperature  profiles  compare  with  the  numerical  ones.  The  comparison  is  then 
extended  to  the  heat  feedback  functions  $g(Ts,P),  the  burning  rates  M(P) ,  and 
finally  the  temperature  sensitivities  a  (tq) . 

A.  Flame  Structure 


Sample  comparisons  are  made  in  Figs.  1-4  for  a  high  Eg  data  set  (RDX1)  and 
one  with  low  Eg  (NC3) .  In  both  cases  the  conditions  of  pressure  and  initial 
temperature  are  lOOOatm  and  300°K.  It  is  clear  that  all  of  the  models  give  very 
simplistic  representations  of  the  heat  release  (q(x))  profiles  (Figs.  1  and  3). 
In  Fig.  1  the  tendency  of  the  BDP  model  to  underestimate  the  flame  standoff 
distance1  for  high  activation  energy  is  evident.  The  high  activation  energy 
case  (Fig.  2)  causes  the  greatest  divergence  in  predicted  temperature  profiles. 
The  asymptotic  theory  does  quite  well,  however,  at  points  near  the  surface.  As 
expected  the  CHR  and  CTRR  models  give  significantly  better  temperature  profiles 
than  the  flame  sheet  models  for  the  low  activation  energy  case  (Fig.  4). 


B.  Heat  Feedback  vs.  Surface  Temperature 

Figures  5-11  display  the  numerical  and  model  values  of  heat  feedback  (i.e., 
dT 

ax  ^x-+o^  as  a  Eunctlon  oE  surface  temperature  for  the  model  propellants.  At 

sufficiently  low  Ts  all  of  the  models  predict  <J>g  t0  approach  the  maximum  value, 
i.e.,  QGMmg-o.  in  this  limit  the  bulk  of  the  energy  release  takes  place  close 
(compared  with  X/MCp)  to  the  surface.  Increasing  Ts  promotes  two  competing 
effects.  The  gas  reaction  is  accelerated  tending  to  increase  4>g »  but  M  also 
increases  by  the  pyrolysis  law  tending  to  blow  the  "flame"  further  away  from 
the  surface  and  reducing  the  effectiveness  of  energy  release  at  any  given 
distance  from  the  surface  (see  the  discussion  of  Eq.  (6)  in  Ref.  1] .  Eventu¬ 
ally  the  latter  processes  dominate  and  decreases  with  Tg. 

Also  shown  in  these  figures  is  <t>_(Tg)  comPute<1  bX  Eq.  (3)  Eor  To=300°K. 

The  intersection  of  <l>s(Ts»'ro)  and  $g(Ts>P)  occurs  at  the  value  of  Ts  peTtinent 
to  the  conditions  (T0,P).  With  this  solution  value  of  Tg  one  can  compute  the 
burning  rate  via  Eq.(l).  s 


An  important  feature  to  notice  in  these  figures  is  that  the  CHR  model  gives 
the  correct  value  of  $g  in  the  high  Ts  limit.  (The  CHR  and  QCHR  models  are 
identical  in  this  limit).  Although  the  solution  to  $s=<I>g  can  always  be  forced 
to  occur  in  this  region  by  considering  a  sufficiently  high  value  of  TQ,  the 


fractional  heat  feedback  ( 


vs0 


)  may  be  so  small  there  that  the  gas  phase 


processes  no  longer  affect  the  burning  rate. 


Another  interesting  feature  is  the  way  that  the  BDP  model  exhibits 
similarities  to  both  the  asymptotic  (A)  and  CTRR  solutions.  This  results 
from  the  BDP  model  using  the  delta  function  (high  Eq)  formalism  but  with  a 
flame  standoff  which  is  appropriate  to  the  low  E^  limit.  (See  discussion  in 
Section  I I I B  of  Ref.  1). 
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Fig.  7.  Heat  feedback  vs.  surface  temperature  for  RDX1 
at  1  atm. 
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In  general  the  CTRR  approximation  is  the  poorest  and  is  most  notably  so  in 
the  high  T  region.  In  this  portion  of  the  curve  the  neglect  of  diffusion  has 

P 

almost  no  effect  on  the  numerical  value  of  $  since  is  small  there.  (See 
Section  IIIA  of  Ref.  1).  b  M 

Finally,  Fig.  11  gives  good  evidence  that  the  asymptotic  theory  is  numeri¬ 
cally  accurate  in  the  high  Eq  limit.  The  approach  to  this  limit  is  also  evident 
in  Figs.  7  and  8.  The  asymptotic  theory,  in  fact,  does  quite  well  generally 
although  Fig.  9  and  high  Tq  cases  in  Figs.  S-8  are  exceptions. 

In  short  it  is  clear  that  no  single  combustion  model  that  we  have  considered 

can  be  relied  upon  for  even  approximate  results  under  all  conditions  of  T  and 
P  for  any  set  of  propellant  data. 

C.  Burning  Rate 

The  burning  rate  as  a  function  of  initial  temperature  for  the  NC2  data  set 
is  shown  in  Fig.  12  at  a  pressure  of  136  atm.  The  NC2  data  set  is  obtained^ 

by  a  fit  of  the  BDP  model  to  double-base  propellant  burning  rates.  However  it 

is  clear  that  the  BDP  model  is  particularly  inaccurate  for  the  parameters  ob¬ 
tained  by  this  fit,  overestimating  the  numerical  rate  by  an  order  of  magnitude. 
Note  that  only  the  asymptotic  theory  gives  good  accuracy  for  this  flame-sheet 
case. 

Figures  13,  14,  and  15  indicate  that,  for  the  most  part,  the  pressure- 
index  (exponent)  is  slightly  less  than  v/2  and  that  the  model  algorithms 
approximate  this  value  fairly  well.  Eq.  (61)  of  Reference  1,  which  connects 
the  pressure  index  with  Es  and  Eg,  is  found  to  be  inaccurate  for  values  other 
than  v/2.  The  absolute  accuracy  of  M  depends  rather  strongly  on  the  parameters 
of  the  data  set.  In  general  the  asymptotic  theory  does  best  at  the  highest  EG 
(RDX1)  and  the  CHR  model  does  best  at  the  lowest  Eg(NC3).  As  expected  from 
the  discussion  of  Figures  5-11,  the  BDP  model  leads  to  an  overestimate  of  the 
burning  rate  for  high  Eg-  All  of  the  burning  rates  in  Figs.  8- 1C  are  computed 
for  an  initial  propellant  temperature  of  300°K.  Ke  note  also  that  the  burning 
rates  are  computed  only  at  1, 1C, 100,1000  atm.  For  clarity,  straight  lines  are 
used  to  connect  these  points. 


P’ 


is  defined  by 


C) 


D.  Temperature  Sensitivity  vs.  Initial  Temperature 

The  temperature  sensitivity  at  constant  pressure,  0 
_  ,d£nM.. 

°P  ° 

op 

Like  the  pressure  sensitivity  (=  )  its  usefulness  lies  in  the  fact  that 

it  is  nearly  constant  over  variable  ranges  of  practical  concern.  Though  nearly 
constant, o  has  nonetheless  been  observed^  to  increase  generally  with  increases 
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Fig.  14.  Mass  burning  rate  vs  pressure  for  RDX1 
Initially  at  300°K. 
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in  T  .  Applying  Eq.  (7)  to  the  asymptotic  theory  (Eq.  (6))  one  obtains 
(cf.  Coates11  Eq.  23) 

Cp  =  (1+  2  +  2RT7  ^  Tf  ’ 


which  predicts  °  to  be,a  monotonic  decreasing  function  of  TQ  (since  Tf  =  T0  + 
Qj/Cp).  The  QCHR  model1,  on  the  other  hand,  predicts  op  to  increase  monotoni- 
cally  with  T  .  Figures  16-21  show  the  resolution  of  this  conflict.  In  all  of 
the  figures  the  abrupt  transition  of  the  asymptotic  theory  curve  (A)  to  the  solid- 

E 


phase  control  limit  (given  by 


R(VVcpr 


■)  does  not  derive  from  Eq.  (8),  but 


rather  from  the  condition  that  the  surface  temperature  cannot  drop  below  T  = 
VQs/Cp-  It  can  also  be  seen  that  in  the  CHR  model  cp  does  decrease  withs 
increasing  T0  at  sufficiently  low  T0.  This  behavior  contradicts  the  QCHR 
formula  derived  in  Ref.  1  even  though  most  of  the  conditions  for  the  derivation 
are  satisfied.  The  discrepancy  can  be  traced  to  the  fact  that  the  QCHR  formula 
was  based  on  a  neglect  of  diffusion. 

Figures  18  and  19  for  RDX1  show  a  slight  hump  in  the  numerical  a  (T  ) 
curves  just  to  the  left  of  the  minimum.  The  numerical  integrations  wire0 
recomputed  in  this  region  with  increased  accuracy  for  the  1  atm  case.  The  hump 
remained,  indicating  that  it  is  apparently  not  an  artifact  of  the  numerical 
method. 

Although  the  abscissa  in  Figs.  12-17  is  given  as  T0,  it  could  also  be 
(To+Qs/Cp)  with  equal  validity  since  Qs  is  fixed  for  each  of  the  data  sets. 

However  if  one  assumes  Tc  to  be  fixed  and  increases  in  the  abscissa  to  be  due 
to  higher  Qs,  one  sees  that  the  consequences  for  o_  depend  on  the  data  set.  For 
example  increasing  Qs  by  50  cal/gm  would  decrease  6p  over  the  practical  range 
for  the  RDX1  data  set  but  increase  op  for  the  NC3  case  (at  1000  atm) .  Thus  it 
would  seem  that  if  combustion  models  were  looked  to  for  guidance  in  making 
formulation  changes  to  affect  op,  a  determination  of  the  effective  kinetic 
parameters  for  the  unmodified  propellant  is  first  required. 

IV.  CONCLUSIONS 

The  objective  of  this  study  was  to  assess  the  numerical  accuracy  of  the 
computational  algorithms  employed  by  several  combustion  models  by  direct 
comparision  with  accurate  numerical  solutions  of  the  conservation  equations. 

This  comparison  was  made  for  a  number  of  data  sets  with  widely  varying  kinetics 
parameters.  None  of  the  algorithms  were  found  to  give  reliable  numerical  accuracy 
for  all  data  sets  and  conditions  of  pressure  and  temperature.  Such  inaccuracies 
can  confuse  efforts  to  associate  idealized  combustion  mechanisms  with  real 
propellants  and  give  erroneous  guidance  for  performance  tailoring  attempts. 


Coates,  R.L.,  "An  Analysis  of  a  Simplified  Laminar  Flame  Theory  for  Solid 
Propellant  Combustion ,  '  Comb.  Soi.  Tech. ,  Vol.  4 ,  pp.  1-8 ,  1971. 
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~  Fig.  16.  Temperature  sensitivity  vs.  Initial  temperature 

for  NCI  at  1  atm. 
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In  spite  of  the  poor  reliability  of  these  algorithms  in  general,  some  of 
the  models  were  shown  to  give  accurate  results  under  certain  limiting  conditions. 
However,  quantitative  criteria  for  identifying  these  limits  may  be  difficult 
to  establish  for  a  priori  use. 

Finally,  we  have  reported  numerical  calculations  of  gas-phase  heat  feedback 
for  a  wide  range  of  kinetics  parameters.  It  is  hoped  that  future  algorithm 
developers  will  find  in  them  a  ready  means  for  validating  their  efforts. 
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LIST  OF  SYMBOLS 


A  -  propellant  molecule  label 
Ag  -  gas-phase  reaction  frequency  factor 
B  -  gas -phase  reactant  label 
C  -  gas-phase  product  label 

C  -  specific  heat  for  solid  and  gas  phases 

P 

D  -  diffusion  coefficient 

-  activation  energy  for  gas-phase  reaction 
Eg  -  activation  energy  for  solid-phase  reaction 
M  -  mass  regression  rate  (mass  flux) 

Mq  -  constant  in  pyrolysis  surface  decomposition  mechanism 

mB~°,  mg+°  -  mass  fraction  of  B  evaluated  at  negative  and  positive  sides  of 
the  solid/gas  interface,  respectively. 

N2  -  number  of  moles  of  C  produced  per  mole  of  B  which  reacts 
P  -  total  pressure 

Qs  -  heat  of  reaction  per  unit  mass  for  solid  reaction  (positive  for  exothermic) 
Q_  -  exothermic  gas-phase  reaction  heat  per  unit  mass  of  B 

b 

Q^,  -  total  exothermic  reaction  heat  (solid  ♦  gas)  per  unit  mass  of  A 
q(x)  -  gas-phase  heat  release  per  unit  volume  per  unit  time  at  x 
q  -  value  of  the  heat  release  function  at  the  surface 
R  -  universal  gas  constant 

j?(x)  -  mass  of  B  reacting  per  unit  volume  per  unit  time  at  x 

i?o  -  constant  mass  reaction  rate  in  Constant  Temperature  Reaction  Rate  model 

r  -  linear  regression  rate  of  propellant  surface 

T,T(x)  -  temperature  at  x 

T0,Ts.Tf  -  initial,  surface,  and  flame  temperatures,  respectively 
ft  -  average  molecular  weight  of  mixture  in  gas-phase 
x  -  spatial  variable 


LIST  OF  SYMBOLS  (Cont’d) 

X  -  heat  conductivity  in  gas-phase 

v  -  gas -phase  reaction  order 

p  -  local  mass  density  in  gas-phase 

Op  -  temperature  sensitivity  at  constant  pressure 

-  heat  flux  from  gas  phase  to  surface  supplied  by  gas-phase  heat  release 

4>  -  heat  flux  from  gas  phase  required  by  solid  to  maintain  steady-state 
s  regression 
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Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1.  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
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make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 
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Telephone  Number: 
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